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ABSTRACT 

Wc present spectro-polarimetric observations of several molecular lines obtained with the 
Submillimeter Array (SMA^l toward the carbon rich AGB star IRC+10216. We have detected 
and mapped the linear polarization of the CO 3-2, SiS 19-18 and CS 7-6 lines. The polarization 
arises at a distance of ~ 450 AU from the star and is blueshifted with respect the Stokes /. 
The SiS 19-18 polarization pattern appears to be consistent with a locally radial magnetic field 
configuration. However, the CO 3-2 and CS 7-6 line polarization suggests an overall complex 
magnetic field morphology within the envelope. This work demonstrates the feasibility of us- 
ing spectro-polarimetric observations to carry out tomographic imaging of the magnetic field in 
circumstcUar envelopes. 

Subject headings: stars: AGB and post-AGB; stars: individual (IRC-l-10216, CW Leo); circum- 
stellar matter; polarization; submillimeter: stars 



1. Introduction 

High mass loss during the AGB phase is one of the main contributors to the return of nucleo-synthesized 
material into the interstellar medium. A proper understanding of AGB mass loss is thus crucial for the study 
of the chemical evolution of the Galaxy. However, the exact mechanism or mechanisms responsible for the 
AGB mass loss is still not clear. After the AGB phase, the star evolves towards the Planetary Nebula (PN) 
phase. During this transition, fast winds are launched that interact with the earlier circumstellar envelope 
(CSE) that was created during the AGB. A large fraction of PNe are observed to be aspherical, and the 
origin of the asphericity is attributed to the influence of a binary companion, a disk, a magnetic field, or a 
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combination of these. The exact onset of asphericity is still unknown, and high angular resolution molecular 
line observations indicate that several CSEs of AGB stars already display various degrees of asymmetry. The 
mechanism responsible for the creation of asymmetries is likely closely linked with that driving the mass loss 
and can be directly probed by high angular resolution observations of CSEs. In particular molecular line 
polarization observations are a unique tool to study potential asymmetries in the CSE and/or determine the 
shape of the circumstellar magnetic field. 

The critical role of magnetic fields in star formation has been probed with dire ct observations of polar- 
ized dust contin uum emission toward both low and high-mass star-forming regions (jCirart et allboodbooot 
Rao et al.ll2009l ). Similar techniques are difficult to apply for AGB stars due to the need for extremely high 
angular resolution (< 0."5) and high sensitivity. A large number of studies have already been made of the 
magnetic field induced polarization of maser lines. These studies have revealed magnetic fields are present 
throughout the entire envelope. The SiO, H2O and OH maser observations indicate that magnetic fields 
appear well ordered and the Zeeman splitting indicates the fiel d strengths range from several Gauss close 



to the stellar surface to several mG at a few thousand AU fe.g..lEtpka fc Diamond! l2004t IVlemmings et al 



2OO5I . l201l[ iHerpin et allbood iKemball et Zlbooot lAmiri et al.lboill ). However, masers probe only a limited 
number of lines of sight through the CSE and in most cases it is thus impossible to fully reconstruct the mag- 
netic field morphology throughout the envelope. Furthermore, the most abundant masers are predominantly 
found around oxygen-right (M-type) limiting the a vailable source sample. Non-mas i ng molecular lines are 
however also predicted to be linearly polarized (e.g. iGoldreich. fc KvlafislEisil . Il982l : iMorris et allll985l) to 



the level of a few percent, and can provide more extensive probes of the entire envelope. 

Polarized en iission has however conclusively been detected from for example CO in a nuin ber of star form- 



ing regions (e.g., ICirart et al]|l999t llai et al.|[20oi ICortes et allboosi : 



Forbrich et al 



200a ). The molecular 



line polarization is due to the anisotropic radiation field from the central star imparting angular momentum 
on the molecules or from anisotropic level populations in the molecular magnetic substatcs when coupled 
with a magnetic field. Without a (non-radial) magnetic field, and assuming a spherically symmetric stellar 
wind, linear polarizat ion should be radia l or tangential and should only be detectable at lines of sight away 
from the central star (jMorris et al.lll985l ). Thus, when these criteria are not met, molecular line polarization 
provides a unique diagnostic o f magnetic field morphol ogy and potential non-radial asymmetries in the stellar 
radiation field. Very rccentlv. IVlemmings et al.l (j2012i ) have reported the detection of CO polarized emission 
toward IK Tau. 

IRC+I0216 (CW Leo) is a weU studied AGB star with a high mass-loss rate, 3 x 10"^ Mq yr-\ 
and a t erminal velocity of 15 km s^^ (e.g. lYoung et al.l 119931) . Due to the relatively close distance of 



150 pc IjCrosas fc Mentenlll997r). this star provides an ideal laboratory for studies of circumstellar chemistry 
( Tenenbaum et aLlboid Patel et allboil ). Unlike typical Mira variables, there are no masers associate d with 
this source, except perhaps a transition of SiS which may be a weak maser (IFonfria Exposito et al Ibood . The 
dusty envelope of IRC-l-10216 shows arc-like structures in scattered light ( Mauron fc Hugging 19991 ) which 
extend over more t han 1^ in angular radius. Closer to the star, d ust emission shows asymmetrical structures 



over scales of 2" (jMen'shchikov et al.l I2OOII : iMenut et al.l 120071 ). Optical and NIR (J HK) interferometr 



reveals complex and tirn e-varying structures on subarcsecond scales close to the star (jTuthill et al 



netry 
200d 



2OO5I : IWeigelt et al.ll2002|) . Non-spherical structu res were also revealed i n linear polarization maps produced 
by dust scattering at 0'.'25 resolution in H band (jMurakawa et al.ll2005l) . 



A previous polarization detection in IRC+10216 was reported by (jClenn et al.l 119971 ) for the CS 2-1 
line. In this letter, we present the detection of the linear polarized emission for the CO 3-2 SiS 19-18 and 
CS 7-6 lines in IRC+10216. This is the first time that maps of the linear polarized emission in IRC+10216 
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are presented. 



Observations 



The SMA observations were taken on 2010 November 24 in the compact configuration. The receiver 
was tuned to cover the 330.6-334.5 and 342.6-346.5 GHz frequencies in the lower (LSB) and upper side band, 
(USB) respectively. The phase center of the telescope was RA(J2000.0)= 9''47'"57?38 and DEC(J2000.0)== 
13°16'43'.'70. The correlator provided a spectral resolution of about 0.8 MHz (i.e., 0.7 km s"^ at 345 GHz). 
The gain calibrators were QSOs J0854-f 201 and J1058-f 015 . The bandpass and polarization calibrator was 
3c454.3, which was observed in a parallactic angle range of ^ 120°. The absolute flux scale was determined 
from observations of Titan. T he flux uncertainty wa s estimated to be ~ 20%. The data were reduced using 
the MIRIAD software package (jWright fc Saultlll993l ). The SMA conducts polarimetric observations by cross 
correlating orthogonal circularly polarizations (CP). The CP is produced by inserting quarter wave plates in 
front of the receivers which are inherently linearly po larized. A detail e d des cripti on of the instrumenta tion 
techniques as well as calibration issues is discussed in Marrone fc Rao ( 2008 ) and Marrone et al. ( 20061 ). In 
order to obtain a more accurate polarization calibration, we solve for the leakage solution independently for 
the strongest detected lines (CO and "CO 3-2, H^CN 3-2, CS 7-6 and SiS 19-18) by selecting a frequency 
range of ~ 1.5 GHz centered within 0.1 GHz with respect to the rest frequency of each line. We found 
polarization leakages between 1 and 2% for the USB, while the LSB leakages were between 2 and 4%. These 
leakages were measured to an accuracy of 0.1%. Self-calibration was performed independently for the USB 
and LSB on the continuum emission of IRC-l-10216. All maps were done with Natural weighting in order 
to maximize the sensitivity, which yielded a synthesized beam of 2'.'6 x 1'.'6 with a position angle of PA~ 0° 
(see caption of Figure [2] for more specific values). Significant polarization was only detected in the CO 3-2, 
CS 7-6 and SiS 7-6 lines, so this paper presents and discuss the detection significance for these three lines. 



3. Results and analysis 

Figure [T] shows the Stokes /, Q and U obtained at the positions of maximum polarized intensity for the 
CO 3-2, CS 7-6 and SiS 19-18 lines. Figure [2] shows the polarization maps for the emission of these lines 
averaged over the velocity range that maximizes the polarized emission, which is different for each line. 

The CO 3-2 lines is the strongest line detected in our observations (the line is 50% brighter than the 
other two lines in the shortest baselines of the visibility plane). This li ne emission is known to b e spatially 



very extended, much beyond the primary beam of the SMA antennas (|Truong-Bach et al. 1991). Howeve r 
the SMA filters out the CO emission that arises from structures larger than about 10" jCirart et al.ll2006l ). 
so the detected emission appears relatively compact. Indeed, the CO averaged emission over the blueshifted 
component shown in Fig. [2] has a relatively compact component (with a radius of ~ 2"), surrounded by a 
weak components that extends up 5" from the center at an intensity level of about 5% of the maximum. The 
CS and SiS show also a compact component with similar dimensions, but they lack the weaker and extended 
component. The Stokes / emission of the 3 lines has a similar brightness, ~ 120 K for the CS and SiS and 
150 K for CO. 

The rms noise of the Stokes Q and U emission appears to slightly increase in the channels where the 
line emission is brightest. This effect is seen in the three lines, though at different levels, from a maximum 
increase of 5%, 13% and 20% for the CS, SiS and CO lines, respectively. This increase is probably produced 
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by the residual leakage (estimated to be of ~ 0.1%. lMarrone fc Radl2008l ) as the total intensity is very strong 



in the ecntral channels. The larger increase for the CO may be due because it is the most extended and 
the brightest line (specially the shortest baselines). Taking into account this increase of noise, there is still 
significant emission in the Stokes Q and U maps of the SiS at the ~ 6-cr level, and in CO and CS lines at 
the ~ 5-(T level. The linear polarization maps were computed by using a 3-cr cutoff, where a is the rms noise 
in the map where the polarization is detected. 

The CO linear polarization arises from both the U and Q components, being relatively bright in the 
later. The CS and SiS lines show mainly polarized emission in the U component. Interestingly, the polarized 
emission in the three lines appears to be blueshifted with respect to the total intensity. 

In order to derive the polarization pattern in the plane of the sky, we have computed polarization maps 
with the emission averaged over the velocity range where the polarization intensity is detected (see Figure [2|). 
The polarization degree at the position where it is strongest is of ~ 2% for the CO and SiS lines, and of ~ 4% 
for the CS. The CO 3-2 polarization arises from two spots, one at the center of IRC+10216 and the other 
located ~ 3" to the East. The polarization vectors are oriented roughly North-South, changing slightly from 
a position angle of PA = —11° at the eastern spot, to —25° at the central spot. The SiS 19-18 polarized 
emission arises mainly from the north-eastern quadrant of the IRC-l-10216 envelope (the polarization peak 
is located ~ 2'.'6 from the center's envelope). There are two other small spots, with a polarized emission 
too marginal to be further considered here. The polarization of the main component has a mean position 
angle of about —41°, but the pattern of the polarization vectors appear to form an are, following roughly 
the contours of the Stokes I emission. The CS 7-6 polarization arises from the envelope's south-western 
quadrant (the polarization peak is located ~ 2'.'9 from the center's envelope). The polarization PA pattern 
is quite uniform with an averaged value of ~ 48°. 

The SiS polarization vectors' pattern suggest a radial distribution. Therefore, we have compared the 
polarization vector direction with the expected radial direction (with respect to the envelope center) at the 
position the vectors. Figure |3] shows the difference between the polarization vectors and the radial directions. 
On one hand, the SiS polarization vectors are all almost perpendicular to the radial direction, i.e., they form 
a nearly perfect concentric arc-like pattern with respect to the envelope's center. On the other hand and 
despite the low polarization statistics, this is not the case for the CO and CS polarization vectors. 



Discussion and conclusions 



The detailed analysis of population of the magnetic sub-levels in rotational lines show that the highest 
polarized emission is expected for volume densities similar to the critical density of the observed transition, 
and depending on the transition and on the molecule, the polarization can still be significant even at densities 
ten times higher (jPeguchi &: Watsonlll984l ). This suggests that the polarization detected in the CO 3-2 line 



should arise at volume densities of ~ 10 em so at the outer regions of the shell, whereas the the SiS 19-18 
and CS 7-6 polarization is expected to trace inner regions, at densities of ^ 10^ cm~'^. 

One of the interesting features is that in the three lines the linear polarization is blueshifted with 
respect to the total emission (this effect is more significant in the SiS line). Considering that the envelope 
is expanding, this suggest that the polarized emission is being detected at the side of a shell facing us and 
with the aforementioned volume densities. In addition, most of the polarized emission arises about 3" offset 
(~ 450 AU in projection) of the envelope's center. Thus the optical depth i s probably playin g an important 
role. Indeed, subarcsccond resolution maps in the IR ( Menut et al. 2007; Leao et al. 2006h and HCN 3-2 



- 5 - 



emission Jshinnaga et alJboogl) show that the molecular distribution is asymmetrical. This suggests the 
anisotropy in the radiation field to be a cause for the polarization pattern to be not distributed spherically. 
This is also in agreement with the single-dish detection of the CS 2-1 line po larization towards the center of 
IRC+10216, which suggests that there is a non-radial polarization pattern (jGlenn et al.lll997l) . 



In circumstellar envelopes and for mm molecular lines the Zeeman splitting is much larger than the 
collision and spontaneous rates, even for magnetic fields strengt h of only few ^tG. Therefore, the polarization 
should be aligned parallel or perpendicular to the mag netic field |Kvlafislll98ilMorris et allll985l ). An overall 
radial magnetic field is expected if it is weak enough to be energetically irrelevant, i.e., the magnetic energy 
is significantly smaller than the mechanical energy of the stellar wind: B <C {mvty^^{Dn)^^, where m is the 



mass loss rate, vt the terminal velocity, D the distance of IRC-l-10216 and fl the angular radius (j Glenn et al 



19971 ). Using the measured values in IRC-l-10216 (see § [T]) and at the distance where the polarization is 
detected (~ 3", 450 AU), this condition is satisfied if i? ^ 8 mG. CN Zeer aan splitting obser vations in 



IRC+10216 indicate a strength of ~ 9 mC at a larger distance (2500 AU, iHerpin et al.l |200£ ) . F or a 



solar-type and toroidal magnetic field configuration (B oc r ^ and r ^ , respectively; IVlemmings 



20111 ) the 



expected strength where the polarization is detected would be in the 50 to 300 mG range. Therefore, the 
magnetic field is strong enough to not be radially shaped by the wind. 

The measured polarized vectors can b e either parallel or perpendicular to the projected magnetic field 
direction in plane-of-the-sky (|Kvlafij l 1983f ) . A proper radiative transfer analysis and a model of the physical 
conditions of the envelope (including the magnetic field configuration) is needed to solve for this degeneracy. 
This applies for the measured CO and CS polarization maps. Nevertheless, in the region where the SiS po- 
larization is detected, the polarization pattern is indicative of a radial magnetic field (see Fig. [3]) . Theoretical 
studies show that in an expanding circumstellar envelope with a radial magnetic field, the polarization will 
be parallel (or perpen dicular) to the radial direction if it arises a t radii lower (or higher) than at certain 
impact parameter, Rj (jPeguchi fc Watsonlll984tlMorris et al ] ll985l) . Rj is the radius where the spontaneous 
emission rate for the J — >■ J — 1 transition is equal to the IR absorption rate. According to iMorris et al 



(|1985[ ). in IRC-H10216 the value of Rj for the SiS 2-1 line is i^a 5 x 10^^ cm. The spontaneous emission 
rate increases with J as J^/(2J 4- 1), and the IR absorption rate goes as R~^, so Rig ~ 2 x 10^^ cm. The 
region where SiS 19-18 polarization is detected arises at a radius of ~ 5 x 10^^ cm from the star. Therefore 
the SiS 19-18 polarization pattern is in agreement with the theoretical predictions if the magnetic field is 
radial in the region where SiS the polarization is detected. 

In summary, we have detected and mapped the polarization pattern for the first time in IRC-l-10216, 
through spectro-polarimetrie observations of the CO 3-2, SiS 19-18 and CS 7-6 lines. Although, the data 
obtained so far lack the sensitivity to allow us to make specific predictions on the magnetic configuration of 
the IRC-l-10216 envelope, the polarization pattern measured discards that the magnetic field configuration 
has a global radial pattern (this is only observed locally where the SiS 19-18 polarization is detected), but 
it possibly has a rather complex magnetic field morphology. In addition, the polarization detection in three 
different molecular lines implies that with the higher sensitivity and angular resolution that the Atacama 
Large Millimeter Array is going to provide, it would be possible to carry out spectro-polarimetrie observations 
for a tomographic imaging of the magnetic field in circumstellar envelopes. However, the ambiguity between 
the polarization direction with respect to the magnetic field direction, implies that in order to properly relate 
the polarization pattern with the magnetic field, a complete radiative transfer analysis should be made. 
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Fig. 1. — Spectra of the Stokes / (top, black line), U (center, blue line) and Q (bottom, red line) emission 
of the CO J=3-2 (left panel), SiS J=19-18 (central panel) and CS J=7-6 lines (right panel). For each line 
this spectrum was taken at the position where the maximum polarized emission is detected, after convolving 
the maps with a Gaussian having a FWHM of 4" x 3". 



EA offset (arcsec: J30DO) 



RA offset (arcsec; J3000) 



EA offset (arcsec: J3000) 



Fig. 2. — Color image of the linearly polarized intensity of the CO 3-2 (left panel), SiS 19-18 (central panel) 
and CS 7-6 lines (right panel), overlapped with the contour maps of the / emission for the respective lines. 
The orange bars represent the polarization vectors. The CS and CO maps show the emission at the wlsr 
velocity of —29 km s~^ averaged over 16 kms~^. The SiS map shows the emission at ulsr =-31.5 km s~^ 
averaged over 20 km s^^. The contour levels are 5, 10, 20, 30, 40, 50, 60, 70, 80, 90, 95% of the peak 
intensity. The wedge shows the polarized intensity scale in units of Jy beam"^ . The synthesized beam is 
shown in the bottom left corner of each panel. 
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Fig. 3. — Distribution of the difference between tlie position angle of the SiS, CS and CO polarization 
vectors and the radial direction with respect to the center of the envelope. We have used a Nyquist sample 
of the polarization vectors to compute this difference, excluding the vector closer than ~ 1" from the envelope 
center. 



